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Abstract. Katz and Pavlovic recently proposed a dyadic model of the Eu- 
ler equations for which they proved finite time blow-up in the H 3 / 2+l Sobolev 
norm. It is shown that their model can be reduced to the dyadic inviscid Burg- 
ers equation where nonlinear interactions are restricted to dyadic wavenum- 
bers. The inviscid Burgers equation exhibits finite time blow-up in H a , for 
a > 1/2, but its dyadic restriction is even more singular, exhibiting blow-up 
for any a > 0. Friedlander and Pavlovic developed a closely related model 
for which they also prove finite time blow-up in Some inconsistent 

assumptions in the construction of their model are outlined. Finite time blow- 
up in the H a norm, with a > 0, is proven for a class of models that includes 
all those models. An alternative shell model of the Navier-Stokes equations is 
discussed. 



1. Introduction 

The Euler equations of compressible gas dynamics typically develop singularities 
in finite time from smooth initial conditions. The singularities are known as shocks 
in that context. The simplest example of this phenomenon is provided by the 
inviscid Burgers, or traffic flow, equation 

,„ . du du 

where u = u{x, t) and x are real, t > 0. Let Q{x, t) = du/ dx, then from Q 

(2) dt +u di = ~^ 

By the method of characteristics, we obtain 

< 3 > §-<•■ 

for C = Q(X(t),t), along the characteristic lines x = X(t) such that dX/dt — 
u(X(t),t) with du/dt — for u — u(X(t),t). Equation J3J exhibits blow-up in 
finite time for negative initial conditions. Its solution is 

(4) C - C ° 



l + Cot 

which is singular at t — — l/£o > if Co < 0. 

This singularity leads to finite time blow-up in the H a Sobolev norm for a > 1/2. 
The H a norm of u, denoted is defined in terms of the Fourier transform 
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u(k, t) of u(x, t) as 

(5) ||u|| 2 ffQ = / (l + \k\ 2a )\u(k,t)\ 2 dk. 

Jr 

Now, the method of characteristics provides the solution of with u(x, 0) = f(x) 
in implicit parametric form as u = f(rj), x = r\ + tf(rj), thus 

2iru{k,t) = J u(x,t)e~ lkx dx 

(6) = / f( v )e-^+tnv)] [l+tf'(T,)]dri 



•ifcfo+t/fa)] 



after integration by parts, where firf) = df/dr/. The singularity occurs when 
1 + tf'(rjo) = for some 7/0 • For simplicity, assume that there is only one such 
770 with f'(r]o) < 0, /"(r/o) = and f"'(r)o) = fb" ^ 0, corresponding to the most 
negative initial slope. The asymptotic behavior of u(k,t) as k — > 00, at the time of 
singularity, then follows by the method of stationary phase as 

u( k ,t)~ -i T /'(r/o)e- 4fc[ " 0+t/(l ' o)1 / e -' ifc */o"s 3 /3! ds 

( 7 > ..,„, ,,,,/. 



where Ai(0) = (2tt) 1 J R e ly3 ^ 3 dy is the Airy function of the first kind evaluated 
at the origin. For example, if f(rj) = —rje - ^ 2 then u(k, 1) ~ i 3~ 1 / 3 Ai(0) fc -4 / 3 . 
Since |tt(fc,t)| = 0(fc~ 4 / 3 ) as fc — * 00, it follows that the H a norm JSJ will diverge 
if a > 5/6. For non-generic initial conditions such that /q" = 0, finite time blow-up 
can occur for smaller a. Smooth initial conditions where f(rj) = —77 in a finite 
neighborhood of the origin, for instance, lead to finite time blow-up for a > 1/2, 
since \u(k, 1)| = 0(fc _1 ) as k — > 00 in such cases. 

It is an open question whether singularities can develop in finite-time for the 
Euler equations when the flow is incompressible. This fundamental question is 
related to the phenomenon of turbulence in fluid flows since classical turbulence 
phenomenology, as well as experimental data, suggest that the energy dissipation 
rate in turbulent flows tends to a strictly positive constant as the fluid viscosity 
tends to zero [5]- 

Katz and Pavlovic |Hj and Friedlander and Pavlovic 0] recently proposed dyadic 
models of the incompressible Euler equations and proved that their models exhibit 
finite-time blow-up in the fl" 3 / 2+c norm. Their proof is restricted to0<e<2/3. 
We show that the Katz-Pavlovic model is in fact equivalent to a dyadic model 
of the inviscid Burgers equation. The latter exhibits finite-time blow-up in H a , 
for all a > 0, much stronger even than the unrestricted inviscid Burgers equation 
which shows blow-up for a > 1/2. Some inconsistencies in the development of the 
Friedlander and Pavlovic model are outlined. It is shown that their model and the 
dyadic inviscid Burgers model belong to a class of models that exhibit finite-time 
blow-up in H a , for all a > 0, even for weak nonlinear couplings. 
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2. Dyadic Inviscid Burgers equation 

After the rescaling u — > 2u, the inviscid Burgers equation Q can be written in 
conservative form as 

, , du d(u 2 ) 

Consider periodic solutions in — 2~ J0 7r < x < 2~- J0 7r, for some jo € Z. Let ui(t) be 
the Z-th Fourier coefficient of 



(9) u(x,t)= Mt)e 2klX 

l — — CO 

where the wavenumber ki = 2 J0 l. The Fourier transform of equation (jHJ reads 

diii 
dt 



(10) ^ = -«fc; ^ U n Ul-n, 



n— — oc 



where the time dependency of ui — u/(t) has been kept implicit. Consider odd initial 
conditions so the solution remains odd for all times, u{x,t) = —u{—x,t), \ft > 0. 
Since u(x,t) is real, this anti-symmetry implies that ui(t) is pure imaginary, hence, 
let itiit) = ivi(t) where vi(t) is real and odd in I, so vi(t) = —v—i(t), and we need 
only consider I > 1. Equation l|10|) becomes 

d x 

(11) -TT = k l X! V n v l-n- 



n— — oo 



Now, artificially restrict to I — 2 m , with m > 0, corresponding to wavenumber 
ki = 2i°l = 2 J0+m . Let j — j + m and define cij(t) = vi(t). Using v n = — u_ n , 
equation 1)11(1 becomes 

da,- 

(12) 



dcij 



3. Reduction of the Katz-Pavlovic model 

The Katz-Pavlovic (KP) model is based, formally, on a wavelet expansion of a 
scalar function u(x, t), with x € M 3 , over the set of dyadic cubes in M 3 . This is the 
set of all cubes having sidelength 2~^ with corners on the lattice 2 _:, Z 3 . If Q is 
a dyadic cube of sides 2~ J , then Q is the unique parent cube of sidelength 2~ J+1 
containing Q and C 1 (Q) is the set of all 2 3 children of Q, each having sidelength 
2~ J ~ 1 . For m > 1, C m (Q) denotes the set of mth order grandchildren of Q, i.e. the 
set of all cubes of sidelength 2~ : >~ m that are contained in Q. 

The KP equations are not derived from the Euler equations, they are chosen to 
mimic the energy conserving quadratic nonlinearity of the Euler equations. The 
nonlinear interactions are restricted to local interactions in wavelet space and de- 
signed to push energy to smaller scales. The KP model equations for the amplitude 
uq{£) of the wavelet localized at cube Q of sidelength 2~ J is |B1 eqn. (2.5)] 

(is) ^ = 2 sm u 2 } _ 2 w)/* UQ £ UQI . 



4 



FABIAN WALEFFE 



The factors 2 5 -?/ 2 were chosen based on the scaling properties of the Euler nonlin- 
earity in R 3 (see [B] for details). 

The KP model is a tree, with each mode having 8 children, however it is a simple 
tree where each edge at the same level has the same weight and each cube interacts 
only with its unique parent and its own 8 children. Therefore if, for any cube Qo 
of sidelength 2~ J0 , the initial conditions are such that all mth order grandchidren 
of Qo have the same amplitude, i.e. if UQ(to) = v,j(to) for all cubes Q of sidelength 
2~ J in C m (Qo) (so j = jo + m), for all m > 1, then they remain so for all t > to- For 
such initial conditions, the dynamics of the Qo branch reduces to a chain model 
duj (t) 
~~dt 

for all j > jo. The factor 2 3 in the 2nd term on the right-hand side arises from the 
fact that each dyadic cube has 2 3 children Q' in C 1 (Q) and we now have UQ>(t) = 
Uj+x(t) for all dyadic cube Q' with sidelength 2~ J . Now if the initial conditions 
are also zero for all cubes with j < j except for the cube Qo, then the complete 
dynamics is described by the chain l|14|) for j > jo together with 

(15) = -2*2^+V/\ u j0+1 , 

for the amplitude of the root cube Qo- 

For such initial conditions, the total energy reads 



(14) ^fZ = 2^"/2 u 2_ i _ 2»2 6 «+ 1 )/ a u J -« i+ i J 



(16) E = J2 2 3m ^ 2 o+m = 2-3^0 £ 2 3 ^ 

m=0 j=j 

The total energy is conserved by the dynamics l|14|) . I|15|l . Our interest is in finite 
energy solutions. 

Define a-j = 2 3 ^ 2 Uj then the energy reads E = 2 -3 -' Y^jLj a j an( ^ the equations 
(HU, ijTB). become 

1 dcij 



— ~ 2 J a, 2 _i - 2 j+1 aja j+1 , j > j , 



30 = -2^ +1 a n 



8 dt 

These equations are identical to the dyadic inviscid Burgers equations l|12f> . after a 
rescaling of time to eliminate the 1/8 factor on the left-hand side. 



(18) 2-^- = 2^l 2 a)_ x - 2^ +1 ^ 2 a 3 a J+1 . 



4. Friedlander-Pavlovic Model 

The Friedlander-Pavlovic model consists of the chain of ODEs 01 e Q n - (3-10)] 
.daj 
~dt 

This scalar chain model is deduced from a more complex vector model. The latter 
vector model is constructed by drawing elements from the Katz-Pavlovic model, 
the Dinaburg-Sinai model |2] and the general philosophy of shell models designed 
to study some features of homogeneous turbulent flows (see for a recent review 
of shell models) . 

The Dinaburg-Sinai model has been discussed elsewhere |Hj- Briefly, it is a model 
of the incompressible Navier-Stokes equations in an unbounded domain, deduced 
by assuming that the nonlinearity is dominated by highly non-local interactions 
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in Fourier space. In that model, small scales do not interact directly with each 
other, they evolve due to the distortion by an infinitely large scale background 
flow with uniform gradient. Strictly speaking, this leads to a linear model of the 
Navier-Stokes equations since the large scale flow would have to be determined 
externally. However, Dinaburg and Sinai inconsistently define the infinitely large 
scale gradient as the linear superposition of the small scale gradients at the origin. 
This leads them to a quasilinear system of equations for the small scales (see |B] for 
details). The derivation of the model is physically and mathematically inconsistent, 
leading in particular to a lack of conservation of energy since the infinitely large 
scale flow contains an infinite amount of energy, but it is possible to prove that 
the resulting system of equations has smooth solutions for limited classes of initial 
conditions 0, 0- 

Friedlander and Pavlovic begin with the shell model approach of considering 
an infinite sequence of nested shells Sj in Fourier space with exponential spacing: 
Sj = {k £ M 3 : 2 J ~ 1 < \k\ < 2 J+1 }, where j is an integer. They consider a 
single time-dependent wavevector kj(t) £ R 3 with associated velocity Vj(t) £ K 3 
in each shell Sj and define the nonlinear interactions following the Katz-Pavlovic 
local interactions model but with some of the structure of the Dinaburg- Sinai non- 
local interactions model. As a result, their wavevector evolution corresponds to the 
distortion by an infinitely large scale flow with uniform gradient, but the latter is 
now defined as the symmetric part of the gradient due to the next smaller scale at 
the origin ([SJ eqn. (2.35)]). This is still physically and mathematically inconsistent 
from a modeling and asymptotics standpoint. It is unclear here why they choose 
to symmetrize their matrix B T which should represent the "large scale" gradient. 
However, since they assume that each wavevector is only distorted by the next 
higher wavenumber, the wavevector dynamics is relatively simple. This allows them 
to construct special classes of initial conditions for which the wavevector dynamics 
and the entire vector structure of the model disappear, except for a dynamically 
insignificant factor of 2 on the left-hand side of JTHJl . The symmetrization of B T 
appears key to that scalar reduction, since if it was defined as a velocity gradient, 



The Friedlander-Pavlovic scalar model H18|l is structurally similar to the dyadic 
inviscid Burgers model \12\. which we showed to be equivalent to the Katz-Pavlovic 
model for a special class of initial conditions. However, the interactions coefficients 
in (|18|) are 2 5j / 2 instead of 2 J . This appears to be another inconsistency of the 
Friedlander-Pavlovic model since they consider a single wavevector of magnitude 
2 J in each shell Sj but their estimate for the nonlinearity, which is valid for wavelets 
in M 3 , is an overestimate for a single Fourier mode. 




would equal a^ +1 ^ /\/2 instead of 



a 



U)/2. 



5. Finite-time blow-up in H' 



Here, we consider the class of chain models 



(19) 




'3-1 



A j+i 



'jo a jo+li 



ajd j+1 , j > j 
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for some integer jo, where A > 1 and the aj's are real. This includes the dyadic 
Burgers model ()12|l when A = 2 and the Friedlander-Pavlovic model l|18(l when 
A = 2 5 / 2 . 

The total energy of system (THfl) is defined as 



(20) E = ^ Ej <ac, Ej = a). 

3=30 

It is easily verified that the total energy is preserved by the dynamics (|19() . thus 
E{t) — E(0), for all t > 0. It is easily verified also that if the initial conditions 
are non-negative, a,j(0) > 0, Vj, then aj(t) > 0, Vj, Vi > 0, since daj/dt > if 
fflj(t) = 0. We consider such non-negative initial conditions hereafter. 

A key characteristic of model (|19|) with non-negative initial data is that the en- 
ergy flux is strictly toward higher wavenumbers (i.e. smaller scales). More precisely, 
let 

OO 

(21) E m =J2Ei 

1=3 

then from lTH?|) 

(22) j t E m = 2XE ] _ 1 a J > 0, 
since a,j(t) > 0. 

The following proof of finite-time blow-up is a generalization of Lemma 5.2 and 
Theorem 5.3 in 0] which are themselves a rewrite of lemma 3.2.2 and corollary 
3.2.3 in 0. 

Lemma 1. For some q with A~ 2 < q < 1, and for some j sufficiently large, if 

(23) E BU) (t )>q j 

— 1/2 

then there is p with (A 2 g) < p < 1, and t in [^^o + P 3 ]} such that 

E B(j+1) (t)>qi+\ 



Proof. The proof proceeds by contradiction. Assume that (|23|l holds but that 

(24) E BU+1) (t) < q> +1 

for all t £ [to, to + pi] and for all p with < p < 1. 

Since E B ^(t) is increasing with E B ^(to) > q 3 , we have 

E m (t) = Ej (t) + E B(j+1) (t) > qi , Vi > t 
and the assumption l|24|) then implies that 

(25) E,(t) > (1 - 9 y , Vi £ [t , t + p>]. 

Now, integrating written for E B ^ +1 ^(t), from t to to + p> and using yields 

/■to+P J 

^Bfa-+i)(*o + (?) = E B(j+1) (t ) + 2A 3+1 / E jaj+ xdt 

(26) 



to 



fto+p 3 

> 2A J+1 (1 - g)<j* / a j+1 dt. 

J t 
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The assumed upper bound flty on E^(j+i) (to + p 3 ) then gives 

fta+P 3 
Jto 



(27) 



a,j+idt < 



q 



2(1 - q)\3+ 1 ■ 



Next, consider the equation for Oj+j.. Since, Ej = a|, it reads 



(28) 



—|— = A J+ ^ - A J+ o i+ io J+2 . 



Integrating from io to to + p> and using the lower bound (|25|) on Ej, the upper 
bound H27fl together with the upper bound <Zj+2 < ^/-Eey+i) < a^ 1 )/ 2 which 
follows from (|24|) . we obtain 



(29) 



a J+ i(io + P J )>g ' +1)/2 



A' 1 /I'l/' 2 KJ- 



gVa 2(1 - q) 



This is in contradiction with the assumed bound aj + i < yEmJ^^ < q 
provided 



(j + l)/2 



(30) 



W 2 > 1, 



for some g and p with < g < 1 and < p < 1, and sufficiently large j. Since 
A > 1, it is always possible to pick q with A~ 2 < q < 1 such that Ag 1 / 2 > 1 and 
then p such that A -1 ;;" 1 / 2 < p < 1, so that condition 1)30(1 can always be satisfied. 

□ 



Define the H a norm of the solution as 

oc 

(si) iic^E^+^n 



]=3a 



where p 3 is the wavenumber associated with amplitude a,j, with /i > 1. For model 
(|19fl . we should pick /i = A, however Friedlander and Pavlovic picked p — 2. A 
proof of finite-time blow-up in H a follows from repeated use of lemma ^ 



Theorem 1. If aj(t) is a solution to system i iy\) with aj(0) > 0, Vj and 
(32) ^ ( ./)(0)>g J 



for some sufficiently large integer J > jo and A 2 < g < 1, then the H a norm of 
the solution becomes unbounded in finite-time, for any a > 0. 

Proof. Applying lemma ^ there is a p with < p < 1 and a time t\ £ [0, p J ] such 
that 

£*?(./+!) (<l) ><Z J+1 - 

Iterating this argument, there is a time tk E [tk-i,tk-i + p J+k ~ 1 ] such that 

E B{.J+k)(tk) > q J+k - 
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Then 

oo 

II 1 1 2 v. 2a/ 2/. \ 

\\a\\ Ha > A* «; (*k) 

l>.J+k 

oo 

(33) > M 2a( J+ *) £ a 2 {tk) 

l>J+k 

>p 2a{J+k) E B[J+k) {t k ) 

> 2a{J+k) J+k 

which blows up as k — > oo provided 

(34) M 2Q g > 1. 
This will occur in finite time since 

(35) t k <p J [l + p + ... + / ,*-i] =p J ^-, 

1 J 1 — p 

which is finite as k — > oo since < p < 1. 

To verify that l|34|) is compatible with the condition 1|30[) required by lemma ^ 
pick q = ^- 2a + s with < 8 < 2a so that q < 1 and (|34l) is satisfied, since > 1. 
Condition l|3UI) requires X 2 q = \ 2 p~ 2a+s > 1 since we need p 2 < 1. Hence, finite 
time blow-up for p — X requires q = p~ 2a+s with 

(36) max{0, 2a - 2} < 8 < 2a, 
while if p = 2 with A = 2 r , r > 0, as in 0], blow-up requires 

(37) max{0, 2a - 2r} < 8 < 2a. 

In either case, finite time blow-up is achievable for any a > 0. 

□ 

This proves that energy-conserving systems of the form i|19|) , with A > 1 , blow-up 
in finite time in the H a norm for any a > 0. This applies to the Friedlander-Pavlovic 
model, for which A = 2 5 / 2 , and is much stronger than the finite-time blow-up in 
jjS/2+e p rove( j m |j Friedlander and Pavlovic concluded that if A = 2 r , then 
r > 3/2 was sufficient for blow-up, but finite time blow-up in H a occurs in fact 
for any r > 0. These results also apply to the Katz-Pavlovic dyadic wavelet model 
which is identical, for special classes of initial conditions, to the dyadic inviscid 
Burgers model (|12f) for which A = 2. Indeed, in the notation of @ and section [31 
with T>j denoting the set of all dyadic cubes of sidelength 2~ J , 

oo 

\\u\\% a >^2^ £ u% 
j=jo QeVj 



(38) > J2 22a3 2 



3Ci-A.) u 2 



J =30 

oo 



2 -3 JO £ 2 2 "',/j 
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and this last sum can diverge in finite time by theorem^ Friedlander and Pavlovic 
choose A = 2 5 / 2 , q = 2~ 3_e and p = 2~ e with e > 0, so q < 1 and p < 1 as required, 
then condition l|3l)fl reads 

(39) X 2 p 2 q = 2 5 - 2e - 3 - £ > 1 

which requires < e < 2/3. Katz and Pavlovic, in effect, choose A = 2, q = 2~ c 
and p = 2 _e , since they work with Uj instead of ctj = 2 3 ^ 2 Uj, leading to the same 
restriction < e < 2/3. 

6. Concluding remarks 

The dyadic models of the three-dimensional Euler equations that have been 
proposed and studied by Katz and Pavlovic [6 and Friedlander and Pavlovic 0] 
are in fact much more closely related to the one-dimensional Burgers equation. 
For the class of initial conditions discussed in section [3] the Katz-Pavlovic model 
is identical to a dyadic inviscid Burgers equation. This is the model obtained 
by restricting nonlinear interactions in the one-dimensional Burgers equation to 
dyadic wavenumbers k — 2 J ko, j integer. The Friedlander-Pavlovic model reduces 
to a similar system but with even stronger nonlinear couplings. The dyadic inviscid 
Burgers equation, and related systems, exhibit finite time blow-up in H a , for any 
a > 0. This is stronger than the finite time blow-up for a > 1/2 that can occur 
in the unrestricted inviscid Burgers equation. Nonlinear steepening is therefore 
enhanced by the restriction to dyadic wavenumbers. 

The idea of studying smoothness and finite time blow-up questions on simplified 
models of the Euler and Navier-Stokes equations has merits. However, one should 
probably consider some of the simplified models that have been studied in the 
turbulence literature (see e.g. £Q). 

An early model, proposed by Obukhov |7j, bears a strong resemblance to the 
dyadic model l|19[) discussed in this paper. The Obukhov model is 

(40) W^^tt = A '"' - AJ+1 «?+i - 

where Vj > is a wavenumber-dependent viscosity coefficient, for instance Vj — 
v\ 2 i , with v > to model Navier-Stokes. The dynamically inconsequential A~ 1//3 
factor on the left-hand side has been introduced so that both models 11190 and l|4(J|) 
have the same power-law steady states 

(41) a 3 = X-^S 1 / 3 A- J '/ 3 

in the unbounded, inviscid limit (y = 0, jo - * — oo), where 

(42) £ = \i E m> 

is the energy flux which is independent of j for the power-law (|41ll . This power law 
is analogous to the Kolmogorov-Obukhov fc~ 5 / 3 power law for the energy spectrum 
in three-dimensional Navier-Stokes turbulence. However, the dynamics of the two 
models are quite different. Model l|19[) is a nonlinear steepening model where energy 
is pushed ever more efficiently to larger wavenumbers. Model l|40|) is an instability 
cascade model. Consider, for example, initial conditions such that all modes are 
initially zero except a;(0) > 0, for some integer I. The nonlinear steepening model 
(|19|l directly starts pushing energy to j > I, but the energy would cascade to j < I 
in the Obukhov model. The energy will cascade to j = l + l only if <z/+i(0) ^ and 
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X l+1 ai > Vi+\. The cascade will continue if ai+2 is non-zero and A i+2 a;+i > or 
in other word if the local Reynolds number is larger than 1 . It would be interesting 
to show what type of finite time blow-up is possible in the inviscid Obukhov model 
and see whether there is any connection between blow-up and the power-law 11411 . 
A connection is expected because a,j > is required for energy cascade to larger j 
and dj(t) is guaranteed to stay positive as long as dj—idj > Aa 2 +1 , when v = 0. 
For power law scaling, <ij oc A™, this requires (3 < —1/3. 
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